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Summary. Effects of stellar rotation on adiabatic oscillation frequencies of ft Cephei 
star are discussed. Methods to evaluate them are briefly described and some of the 
main results for four specific stars are presented. 

i-C ■ 1 Introduction 

O ■ Main sequence (MS) massive stars are usually fast rotators and their fast rotation 

f~i | affects their internal structure as well as their evolution. The issue which is 

adressed here is what information can we obtain- about rotation - from the 
oscillations of these massive, main sequence stars ? 

The following seismic diagnostics for rotation using non axisymmetric modes 
will be discussed: 1) rotational splittings as direct probes of the rotation profile. 
More precisely, we study the effects of cubic order in the rotation rate compared 

£vqj | to effects of a latitudinal dependence of the rotation on the splittings; 2) splitting 

asymmetries as a probe for centrifugal distorsion. The case of 3) axisymmetric 

C*~) ' modes as indirect probes of rotation throughout effect of rotationally induced 

[^. ■ mixing on the structure will also be considered. 

<^2 \ Results discussed here are obtained with perturbation methods. For nonper- 

turbative methods and results, we refer the reader to Lignieres et el. (2006), 
Reese et ai. (2009) and references therein. 

The paper is organized as follows: in Sect. 2, properties of pulsating B stars 
are recalled with emphasis on the uncertainties of their physical description that 
can be addressed by seismic analyses. Sect. 3. recalls the theoretical framework 
for seismic analyses of relevance here. In Sect. 4, seismic analyses of four /3 Cep 
are presented. In Sect. 5 a theoretical study compares the modifications of the 
rotation splittings due either to latitudinal dependence of the rotation rate, f2, 
or to cubic order (0(/2 3 )) frequency corrections. Some conclusions are given in 
Sect.6. 



> 

ON 



* 



Mariejo Goupil 



^ 

■^ 



*< 


i 












"*"■*. 


. 


J 




/ 




y 


*"** ^ 


■^^y' 




~y~l ~^ 


/ 


/ 


^y ^ 


/ 


^y*^ ^ y 


-s ^ 


f*y r'' 1 ' 


/ '">< 


," 


/ s*-' 


"'"■'-. 


f/£' 


_-.--'' 


f/' ,--'■■" 




l& 




F y 




/I / 




fi'' r 


V^IQO km a-i 


ft 


V M =3Vt 


ft 


overshooting 


^s i 


no overshooting 










4 35 4.4 



Fig. 1. Evolutionary tracks for 9Mq models with neither rotation, nor overshoot 
included (dotted line), with overshoot included but not rotation (short dashed line) 
and with rotation included but not overshoot (long dashed and solid lines) (from Talon 
et al., 1997) 



2 Massive main sequence stars 

0-B stars are characterized by a convective core and an envelope which is es- 
sentially radiative apart a thin outer region related to the iron opacity bump. 
Important incertainties regarding the structure and future evolution of these 
stars are: 

-the extent of chemical clement mixing beyond the central instable layers as de- 
fined by the Schwarzschild critcrium 

-Transport of angular momentum because the rotation can play a significant role 
in chemical element mixing 

Convective core overshoot: In ID stellar models, the convective core is delim- 
ited by the radius r zc according to the Schwarzschild criterium V a d = V ra( j. How- 
ever this corresponds to a vanishing buoyancy force: the eddies are then strongly 
slowed down but still retain some velocity. Hence due to inertia, eddies move be- 
yond the Schwarzschild radius till their velocity vanishes that-is over a distance 



d ov such that the effective convective core radius becomes r 



De- 



spite theoretical investigations (Zahn, 1991, Roxburgh, 1992), the overshooting 
distance computed in ID stellar evolutionary models usually remains a rough 
prescription i.e. it is assumed that d ov = a ov mm(r zc , H p ) with H p the local 
pressure scale height and a ov is a free parameter. Empirical determinations from 
observations (Schaller et al. 1992, Cordier et al. 2002; Claret, 2007 and refer- 
ences therein) yield a wide range for a ov , namely [0-0.5] H p . The adopted value 



for this free parameter has important consequences for the evolution of a model 
with a given mass: with a higher luminosity, it is older at given central hydro- 
gen content (X c ) on the MS and reaches the end of the MS with a larger mass 
core- total mass ratio. On a statistical point of view, the value of a ov affects the 
thickness of the MS on a HR diagram as well as the isochrones. Core overshoot 
has therefore an influence on stellar age determination (Lebreton et al. 1995, 
Lcbreton, 2008). 

Rotationally induced mixing in radiative regions: Departure from thermal 
equilibrium generated by the oblateness of a rotating star causes large scale mo- 
tions, the meridional circulation. As differential rotation also induces turbulence, 
competition of these two processes can result in (rotationally induced) diffusion 
of chemical elements (Zahn, 1992 and subsequent works). The evolution of a 
given chemical specie j with concentration Cj is governed by a diffusion equation 
(for a review, Talon, 2008; Decressin et al. 2009): 

dc,j 1 d , 9 , 1 d , , _ dej . , . 

at r z or r z or or 

where the first term represents nuclear transformation and the second term 
atomic diffusion with Vj_ p the diffusion velocity of particles j with respect to 
protons and where the tuirbulent diffusion coefficient D t — D e // + D v , D e ff 
comes from the meridional circulation and D v from the turbulence. As -D e // de- 
pends on the vertical meridional velocity U r , chemical and angular momentum 
evolutions must be solved together. Hence one also solves an (diffusion-advection) 
evolution equation for the angular momentum : 

dr 2 fl 1 d , 4 ^ TT . 1 d . 4 dQ, 

where v v is the vertical turbulent viscosity related to rotational instabilities. The 
current picture is that the vigor of the meridional circulation is controlled by 
the magnitude of the surface losses of angular momentum. Hence for hot, high 
mass stars which lose mass but much less angular momentum, one expects no 
efficient angular momentum internal transport. The rotation profile then essen- 
tially results from expansion and contraction within the star during its evolution: 
i.e. high ratio of core rotation over surface rotation. This is well reproduced by 
rotationally induced mixing of type I (Talon et al. 1997). On the other hand, for 
cool stars with extended convective outer layers, dynamo generates an efficient 
magnetic driven wind which is efficient to drive important angular momentum 
losses and internal transport. This mechanism however is not sufficient enough in 
the solar case to make the observed rigid rotation in the radiative solar interior 
and one must calls for to other mechanisms (waves, magnetic field) (see Talon, 
2007; Rieutord, 2006 for reviews). This shows that many open questions related 
to stellar internal rotation and its gradients subsist. An important issue then is 
to locate regions of uniform rotation and regions of differential rotation (depth 
and/or latitude dependence) inside the star (f2 core / f2 sur f). Another problem 
which must be solved is to disentangle effects of overshooting and rotation on 
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Fig. 2. left: HR diagram and instability strip for beta Cephei stars. Full dots represent 
confirmed /3 Cephei stars and open dots : candidates. The dashed lines delimitate the 
instability strip for the fundamental radial mode (adapted from Stankov & Handler 
2005). 



mixed central regions and extension of convective cores. Indeed the rotationally 
induced chemical mixing affects the evolution of the star, its internal structure 
and oscillation frequencies as does core overshoot although in a different way 
(Talon et al, 1997; Goupil & Talon, 2002; Montalban et al, 2008; Miglio et al, 
2008; Thoul, 2009). FigJT] illustrates the respective effects of element mixing by 
core overshoot and rotation on the evolution of a 9 Mq main sequence model in 
a HR diagram. 

Seismology of O-B stars can bring some light about these processes. More 
specifically, (3 Cephei stars are good candidates for this purpose (Montalban et 
ai., 2008; Miglio et al, 2008; Goupil & Talon, 2008; Miglio et al, 2009; Lovekin 
et al, 2008; Lovekin & Goupil, 2009). Indeed, unlike S Scuti stars, (3 Cephei 
stars do not present significant outer convective layers which makes the mode 
identification more trustworthy provided the star is slowly rotating or that its 
fast rotation is taken into account in the mode identification process (Lignieres 
et al, 2006, Reese et al, 2009; Lovekin et al. 2008; Lovekin et al 2009). 



2.1 (3 Cephei stars 

f3 Cephei stars are main sequence stars with masses roughly larger than 5 — 7M Q 
(Fig. [2]). They oscillate with a few low degree, low radial order modes around the 
fundamental radial mode i.e. with periods around 3-8 h. The modes are excited 
by the kappa mechanism due to the iron bump opacity. These pulsating stars are 
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Fig. 3. Propagation diagram for model A: a 8.5 Mq model with T e // = 22230^ and 
initial X = 0.7, Z = 0.019 (no rotation, no overshoot). The Lamb frequency (dashed 
line) is plotted for £ = 1. Normalized squared frequencies discussed here are found in 
the range er 2 =5-15. 



located at the intersection of the main sequence and their instability strip in a 
HR diagram (Figj2]). For more details about (3 Cephei stars, we refer to reviews 
by Handler (2006), Stankov & Handler (2005), Pigulski (2007), Aerts (2008). 

Sofar the observed modes have been identified as pi, p2, gl modes. We recall 
that p modes are propagative when uj 2 > N 2 and uj 2 > Sf (for more details, 
see Christensen-Dalsgaard, 2003 CD03). The squared Brunt- Valssala (buoyancy) 
frequency is defined as 



N 2 = 9_ ( l dln P 



dlnp 
i"i d In r d In r 



(3) 



with p, p, g respectively the pressure, density and gravity of the stellar medium 
and Pi the adiabatic index. The squared Lamb frequency is defined as 



Sf = (k h c s ) 2 = £(£ + !)- 



(4) 



with kh the horizontal wavenumbcr of the pulsation mode and £ the degree of 
the mode (when its surface distribution is described with a spherical harmonics 
Y™{0, (/>)). The local sound speed c s is given by: 



(?r 



(5) 



For g- modes, the propagative region is delimitated by uJ 2 < N 2 and u 2 < Sf. 

For P Cephei stars, mixed modes propagate as g mode in the inner part and 
as p mode in the outer part of the star. Depending on the evolutionary stage of 
the star, one expects some of the detected modes to be of mixed p and g nature. 
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Fig. 4. Histogram of projected rotational velocities for /3 Cephei stars (from Stankov, 
Handler 2005). 



Modes with frequencies around that of the fundamental radial one (normalized 
frequency a = u/Q K ~ 2 - 3 with Q K = (GM/R 3 ) 1 / 2 , R the radius and M the 
mass of the star) can be mixed modes. This can be seen in FigJ3] which shows a 
propagation diagram for a typical case, model A, a model with a mass 8.5 Mq 
and an age = 19.9 Myr, a solar metallicity to hydrogen ratio Z/X — 0.019 with 
X = 0.7 and log T eff = 4.347 and log L/L Q = 3.723 that therefore lies in the 
middle of the main sequence and instability strip for these stars. 

Rotation of j3 Cephei stars ranges from slow (rotational velocity v < 50 km/s) 
to extremely rapid (v > 250 km/s) (FigHJ). Effects of uniform rotation start to 
modify significantly the tracks in a HR diagram beyond v — 100 km/s for these 
masses (Lovekin et ed., 2009). For v = 100 km/s, with a stellar radius R = 
4.94i? Q , model A is characterized by Q/(2 K ~ 0.175 where Q K = {GM/R 3 ) 1 / 2 
is the break up angular velocity. 



3 Theoretical framework 

In this section, we recall the theoretical framework within which seismic observa- 
tions of these stars can be interpreted in terms of rotation (for more details, the 
reader is referred to Goupil (2009) and references therein). For sake of notation, 
we recall first the non rotating case. 



3.1 No rotation 

Adiabatic pulsation studies consider the linearized conservation equations for a 
compressible, stratified fluid about a static equilibrium stellar model character- 
ized by Po,po,Ti, 4>q respectively pressure, density, adiabatic index, gravitational 
potential profiles. The equation for hydrostatic equilibrium is : 
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Vp = ~Po V0 O (6) 

Assuming the fluid displacement <5r(r, t) of the form <5r(r, t) = £(r) exp(iu>ot), 
the linearized momentum equation then is: 

£ (£)-A)^U = (7) 

with 

where £q is a differential operator acting on £; p',p',<f>' are the Eulerian per- 
turbation for the pressure, density and gravitational potential respectively. One 
must add boundary conditions (Unno et al, 1989) and this gives rise to an eigen- 
value problem where luq is the eigenvalue for the nonrotating case and £ is the 
eigenfunction for the fluid displacement. In the following, we will keep the nota- 
tion: v in /iHz or c/d ; lu in rad/s; a = ^/(GM/R 3 ) 1 ^ 2 the normalized frequency. 
One defines the scalar product: 

< a|b >= f a* • b d 3 r (8) 

Jv 

where * means complex conjugate and where V is the stellar volume. The scalar 
product of £ with Eq[7]then yields: 

< C|£o(€) - Po«& >= f C • (A)(€) - A>wfc) d 3 r = 

The eigenfrequency can be obtained as an integral expression: 

ul = j < Z\jr (S) > (9) 

l = \J r-(Vp' + poV0' + p'V0 o )d 3 r (10) 



or 

with 



/= / (r-^d'r (11) 

Jv 

In absence of rotation, the eigenmode displacement is written in a spherical 
coordinate system with a single harmonics, Y™(9, 4>), with a spherical degree £, 
an azimuthal number m being the number de nodes along the equator 

£(r) - Cr(r) YT e, + fr(r) V h Y e m (12) 

where the first part is the radial component and the second term the horizontal 
component of the fluid displacement. The horizontal divergence is 

„ d Id 



d9 v sm9dcf> 
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The divergence of the fluid displacement is written as: 

V • f = A Y e m (13) 

with 

A = --r^ - - & 14 

r z ^ r r 

and A — £(t + 1). The perturbed density p'(r) = p'(r)Y™ is given by the lin- 
earized continuity equation: 

p'(r) = -V • ( Po £) = -^ r - PoA (15) 

The perturbed gravitational potential <j>'(r) — (fi'^Yf 1 is given by the perturbed 
Poisson equation: 

V 2 <£' = AirGp' (16) 

The pressure perturbation p'(r) is related to the density perturbation p'{r) 

by the adiabatic relation (Unno et aI1989) where S means here a Lagrangean 

variation: 

5p =r Sp 

Po Po 

3.2 Including rotation 

In presence of rotation the centrifugal and Coriolis accelerations come into play. 
The centrifugal force affects the structure of the star - the star is distorted- and 
causes a departure from thermal equilibrium which generates large scale merid- 
ional circulation and chemical mixing. Accordingly, the resonant cavities of the 
modes are modified. The static equilibrium (averaged over horizontal surfaces) 
ID stellar model is modified and characterized by Po,n{r),po,n(r), A, n{r)i 4>a.n{ r ) 
with f2(r, 9) the rotation rate. The Coriolis force enters the equation of motion 
and affects the motion of waves and frequencies of normal modes. The linearized 
equation of motion is modified. As rotation breaks the azimuthal symmetry, it 
lifts the frequency degeneracy: without rotation, 1i-\-\ modes with given n, £, m 
(m = —I, £) have the same frequency ujq (omitting for shortness the subscripts 
n,£). With rotation, the same 2£ + 1 modes have different frequencies w m and 
the rotational splitting is defined as : S m — (u) m — u>o)/(m). One also uses 
S m = uj m — C0m-i an d the generalized rotational splitting: 

S m = ^-^-^ (17) 

2m 

where uj m is the mode frequency. These various definitions are equivalent only 
at first perturbation order in the rotation rate fi; the first two are used when 
only a few components arc available. 



3.3 Rotational splittings 

At first perturbation order in Q, only the Coriolis acceleration plays a role. The 
linearized equation of motion including the effect of Coriolis acceleration (2i7 x v) 
in a frame of inertia is 

A)(£) ~ Po (w m + mQfi + 2p (w m + mQ) Q K.£ = (18) 

with /C£ = ie z x £ and £ is the displacement eigenvector in absence of rotation 
and e z is the vectical unit vector in cylindrical coordinates. The nonrotating case 
is recovered by setting Q = 0. One then expands the displacement eigenfunction 
as £ = £o + £i an d the eigenfrezquency as u> m — Uq—o + wi, m where wn=o, £o 
correspond to the eigenfrequency and eigenfunction for a nonrotating star and 
wi. m , £i give the first order correction due to Coriolis acceleration. Keeping only 
terms up to O(Q), one obtains: 

A)(£i) - Po ^n=o £i ~ 2 Po wn=o(wi, m + mQ) £ 

+ 2 Po un=o n IC^o = (19) 

The correction to the eigenfunction £i can be chosen so that < £o|£i >= 0. 
Taking the scalar product EqJS] of £ with Eq[TjJ] and keeping only terms up to 
O(Q) yields: 

< £o| [ A)(£i) - Po w^ =0 li - 2,00 u n=0 (uj lym + mQ) £ 

+ 2p uj n=Q Q /C£o ] >= (20) 

from which one derives for a mode with given n, £ 

wi, m I = ( &'(fiK- mQ) £ Po d 3 r (21) 

Jv 

which is rewritten as: 

pR pit 

0Jt, m = / K m (r,0) Q(r,0) dOdr (22) 

Jo Jo 

where the analytical expression for the kernels K m is given in Appendix. At first 
order O(Q), the generalized splitting Eg IT7I then is given by 

S m = " 1 ' m :" 1 - m (23) 

2m 

Assuming a shellular rotation Q(r), the splitting S m becomes m independent 
and one has: 

S = f K(r) Q(r) dr (24) 

Jo 

with the rotational kernel 
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K(r) = -j {£ - 2Uh + (A - 1)3) Po r 2 (25) 

and mode inertia (Eg fTTj) : 

/= / (g+MDpo^dr (26) 



S-- 


Jo 


K(x) Q{x) 


dx = 


= fl«A 


+ n e /3 e 


/8c = 


-F 

Jo 


K(x) dx : 


Pe 


-f 


K{x) dx 



with again A = £(£ + 1) and R the stellar radius. For a uniform rotation, this 
further simplifies to 

S = Q ; P = / K(r) dr (27) 

Jo 

f3 is assumed to be known from an appropriate stellar model, S is measured and 
Q is inferred. This will be used in Sect. 4 for j3 Cep stars. 

When only a few measured splittings are available, information about the 
internal rotation is limited so one assumes for instance a uniform rotation for 
the convective core with the angular velocity fi = i? c (for x — r/R < x c ) and 
a uniform rotation for the envelope fi — fi e for x > x c . Both values are the 
unknowns. Inserting into Eql2"4l 



with 



The detection of 2 triplets £ = 1 for instance yields fi c , fi e and fi c / fi e provided 
/3 C , f3 e are given by a model close to the observed star. This type of approach 
was used to determine whether the star is in rigid rotation or not for a 5 Scuti 
star (Goupil et al, 1993); for white dwarfs (Winget et al. 1994; Kawaler et ai., 
1999) and recently for (3 Cephei stars (Sect. 4 below) and SdB stars (Charpinet 
et ai., 2008). 

3.4 Splitting asymmetries: distorsion 

At second order in the rotation rate, the centrifugal acceleration comes into play. 
This has several consequences on the oscillation frequencies (for a review Goupil, 
2009). One is that the split components are no longer equally spaced. It is then 
convenient to define A m the splitting asymmetry as 

A m = Wo - ^(0J m + 0J- m ) (28) 

In order to interpret observed asymmetries, let consider a given multiplet of 
modes (i.e. with specified n,t). Its oscillation frequencies, u m (m = —£, ..,£), are 
computed up to second order 0{Q 2 ) as: 



11 
n 2 2 

w m = cjo.fi + m S\m\ H (Xi + m A 2 ) (29) 

wo.fi 

whereand cJo.fi is the eigenfrequency for a static model including the horizon- 
tally averaged centrifugal acceleration. The second term is the splitting Eq|23l 
due to Coriolis effect and fi is an averaged rotation rate. The last term is the 
asymmetry due to the non spherical part of the centrifugal distorsion which 
dominates for high radial order modes. Expressions for Xi,X 2 can be found in 
Saio (1981), DG92, Soufi et al. (1998); Suarez et al.(2006); Goupil (2009). For 
low radial modes such as those excited in /3 Cep stars, the second order Coriolis 
contributions to X\, X 2 remain significant. According to Eq |29[ the asymmetry 
is then given by: 

A m = (—) m 2 X 2 (30) 

V^o.fi / 

Let consider again the linearized equation of motion including now the centrifu- 
gal acceleration: 

4),fl(*) - Po.fi cj 2 <£ + 2p 0t nwQK{i) + C 2 {£) - p 2 Qj 2 £ = (31) 

where cD = w m +mfi. The spherical part of the centrifugal acceleration is included 
in the spherical ID model, therefore the linear operator depends on the rotation 
rate i.e. 

£ ,fi(O = W + po.fiV^' + p'V^o.fi (32) 

and for the non spherical distorsion 

£ 2 (£) = p'(— Vpo.fi - Vp 2 ) + p 2 V<// + po.fi e s r sin0 V U 2 ■ £ (33) 

Po.fi 

where e s = sin# e r + cos6 eg in a spherical coordinate system (e r ,eg,e0).The 
subscript 2 indicates departure from sphericity p 2 , p 2 , <fi 2 for the pressure, density 
and graviational potential respectively. Again using the scalar product EqjSJ one 
writes 

< £o\£o,n(£) - Po.fi Co 2 £ + 2po.fi u n K(£) > 

+ < ^\(C 2 (£) - p 2 io 2 £) >= (34) 

One then assumes an eigenfunction of the form £ = £o + £i + £2 and the 
eigenfrequency as u) m = cjo.fi + u)\ m + uj 2 where the unknown now is uj 2 . Solv- 
ing Eq[34] for cj2 m leads to an integral expression for X\ , A2 and therefore an 
expression for A m of the form: 

Am = m 2 J i7 2 (x) K 2 {x) dx (35) 

Jo 

where K 2 (x) include effects of distorsion of the structure throughout p 2 , p 2 and 
depend on the eigenfunction £. Detailed expression for K 2 (x) can be found in 
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DG92, Soufi et al. (1998), Karami (2008), Suarez et al. (2006), Goupil (2009). 
An exemple for K 2 (x) is shown in FigJ7]and discussed in Sect. 4. 2 

Splitting asymmetries can provide probes of the internal structure which 
differ from those derived with the splittings S m as the corresponding kernels are 
different. When only a few observed asymmetries are available, one can proceed 
as for the splittings fSect l3.3l above). Assuming a rotation profile of the simplified 
form: 

J? 2 (x) = fi\ for x c < x 

fl 2 (x) = Q\ + 2(x - x c ) f}'£2 c + {x — x c ) 2 Q 2 for x c < x < x e 

Q 2 {x) = n 2 for x e < x (36) 



with 



then 



Q' = Qe Qc (37) 



A-m — 



■{n 2 c p 2 ,o + 2Q'n c p 2 ,i + h,2 n' 2 + fo^t) (38) 

where f2 c , fl' are assumed known from the splittings (Sect. 3. 3) and 

fa, q = I C (x~x c ) q K 2 (x) dx (39) 

h,e = / K 2 (x) dx (40) 



Determination of the f3 2 coefficients then brings some information on the ker- 
nels K 2 (x) with the promising prospect of deriving constrains on the rotationally 
distorted part of the stellar structure. 

3.5 Axisymmetric modes: rotationally induced mixing 

Centrifugal departure from spherical symmetry has important effects on all 
modes including the axisymmetric modes. Indeed the values of the ?ti = mode 
frequencies are shifted when compared to those of non rotating models. Hence 
the differences 






8u = uj - w { Q = 0) = ( - — ) m 2 Xi (41) 



from Eg 1291 between frequencies of a given mode from a model including rota- 
tion and a non rotating model can be an efficient diagnostic for rotation effects 
although some care must be taken when defining the fi — stellar model for 
comparison. This has been extensively discussed in past publications (Chan- 
drasekhar & Lebovitz, 1962; Saio, 1981; Gough & Thompson, 1990; DG92; CD03, 
for a review, see Goupil, 2009). 

Another (indirect) effect of the star oblatness on frequencies, as already 
mentionned in Sect (21 is due to the departure from radiative equilibrium which 
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generates large scale motions (meridional circulation), differential rotation and 
consequently shear turbulence. All this concurs to affect the rotation profile. 
It also causes mixing of chemical elements which affects the prior evolution of 
the observed star and therefore its structure. These structure changes must be 
computed by coupling both evolutions of the angular momentum and the chemi- 
cals, as already mentionned in Sect. 2. These equilibrium structure modifications 
affect all modes as compared to those of a nonrotating star, including the ax- 
isymmetric modes. The effect on the frequencies can be quite significant as was 
discussed by Goupil & Talon (2002) and quantified by Montalban et al. (2008), 
Miglio et al. (2008), Goupil & Talon (2008) (see Sect.4.3 below) 

We consider here only the effect of the structure modifications due to rota- 
tionally induced mixing on the axisymmetric mode frequencies. The Coriolis or 
the centrifugal accelerations then are not included in the linearized equation of 
motion. Hence the linearized equation of motion including rotationally induced 
mixing yields the usual integral expression for the eigenfrequency of a nonro- 
tating model, Eq llOi except for the structure quantities such as the density, the 
pressure, the gravity (p,p 7 g resp. ) etc... which are modified by the rotationnally 
induced mixing. As they now depend on the rotation rate, we write them as 
pn,pn,gn---- The linearized equation of motion including rotationally induced 
mixing in a ID spherically symmetric stellar model then is given by: 



U 0J2 — 

with the mode inertia 



J- I Co' (Vj/ + poV4>' + p'Vfa) d 3 r (42) 

in Jv 

In = I (tn ■ *fl) Pn d 3 r 

Jv 

From now on for sake of shortness, we omit the subscript Q for the eigenfunc- 
tions. We define the dimensionless variables according to Dziembowski (1971) 
(see also Unno et al, 1989): 

m = -; y 2 = — «-' + —); 2/3 = — 0' (43) 

r gnr pa gar 

Starting with Eq l42l integrations over surface angles and a few integrations 
by part for the radial part yield: 

1 f , w __ , __ N d\np n _ ,_ __ ^ 3 



w o,f3 = T~ / (- x (Vi + V2) - -j: 2/1 (2/1 + 2/3)) gn Pn r dr (44) 

In Jr d\nr 

where we have assumed that the surface integrals vanish. From its definition 
EqUl 

A = Vg.n (yi -2/2 + 2/3) 

with 

1 d\np n 



V, 



9,0 



Zio d\nr 
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. Note that there are several alternative equivalent expressions for Wq q . 

Differences between the structure of a model including rotationnaly induced 
mixing and that of a model which does not result in differences in the eigenfre- 
quencies which we note 8u> = wo,j? — wo.j?=o- We will see in Sect. 4. 3 that the 
structures of the models indicate that pa and its derivative with respect to the 
radius, the gravity go, the density pn are not significantly modified compared 
to the derivative of the density with respect to the radius. Accordingly using 
Eq l44l and keeping only the first order terms, one obtains: 



2 cjf2 = o In=o 
I In pn 
dlnr 



/ 5 ( ^ " ) Vi ( yi + y ^ 9n =° pn =° r3dr ( 45 ) 



where we have also assumed that the perturbations of the eigenfunctions yj o — 
Vj,o=o (j = 1,3) are negligible at first order. For massive main sequence stars 
the largest difference 8(d In pa/d In r) arises near the convective core (Sect. 4. 3). 
Largest frequency differences therefore are expected for mixed modes compared 
to p-modes. Note that the same interpretation can be obtained with differences 
in the Brunt- Valssala behavior. Indeed at the same level of approximation, one 
has from Eq(3] 

'(f ) " -»(w) « 

For high frequency (i.e. pure) p-mode which propagate significantly above 
the V/x region, the difference S I d In pnj d In r ) is essentially negative. In addition 
1 2/3 1 << 1 2/i I so that we obtain 

8u> 



2 tun = o Iq= 



J R S{ lm7 ) Vl 9Q =* PQ =° T " dr > ° (4?) 

which is small and positive. For mixed modes having high amplitude in the V/i 
region, 5 1 din pajd In r J can be positive and the frequency difference can be 
large and negative as illustrated in Sect. 4. 3. The difference 8u) is quantified and 
discussed in the case of a f3 Cephei model in Sect. 4. 3. 



4 Seismic analyses of four (3 Cephei stars 

We discuss 4 j3 Cephei stars which have been the subject of seismic analyses 
and for which information about rotation and core overshoot has been inferred: 
V836 Cen (HD 129929); v Eridani, 9 Ophiuchi and 12 Lacertae (see also Thoul, 
2009). Schematic representations of the frequency spectra for the first three stars 
are displayed in FigJS] and FigEl These four stars are relatively slow rotators 
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76 78 80 




Fig. 5. left: Schematic representation of the frequency spectrum of HD129929 (data 
from Aerts et al, 2004)- right: Rotational kernels for the excited pi and gl modes of 
HD129929 in function of the radius r/R normalized to the stellar radius (from Dupret 
et al, 2004). 



(with surface rotational velocities smaller than « 70 km/s). Determination of 
the luminosity, effective temperature and location in the HR diagram for these 
slow rotators are not significantly affected by rotation. 



4.1 Rotational splittings 

HD129929 

is a main sequence ~ 9M Q star for which one £ = 1, p\ triplet has been de- 
tected and identified as well as one radial mode and 2 successive components 
of the £ = 2, gi mode as represented in FigJS] (Aerts et al., 2004, Dupret et 
al., 2004). From the triplet and assuming a solid body rotation, , one uses 
S = Q j3 (Eq [M|) as explained in Sect l3.3l With f3 known from an appropri- 
ate stellar model, the measured splitting for the £ = l,p = 1 triplet S gives 
v r ot = 3.61 km/s but from the two successive components of the £ = 2 multiplet, 
one obtains v ro t = 4.21 km/s, clearly indicating a nonuniform rotation (Dupret 
et al., 2004). Assuming therefore a uniform rotation for the convective core with 
angular velocity fl = fl c and a uniform rotation Q = fi e for the envelope of the 
star, the splittings then obey S = J7 C /3 C + Q e f3 e where /3j are the integral for 
the core or the envelope fSect !3.3[) . It is found that \/3 \ << \/3 e \ that-is actually 
the detected modes do not efficiently probe the convective core. This can be 
seen with the associated rotational kernels in FigJS] which have no amplitude 
in the core. Therefore fl c is taken as the rotation rate of the radiative region 
in the /i-gradient region above the convective core (with fj, the mean molecular 
weight). Assuming a linear depth variation of the angular velocity in the enve- 
lope Q(x) = J?o + {% — xo)n 1: the splittings must obey S = /2oA) + ^i/3i where 
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Fig. 6. Schematic representation of the power spectrum of left: 9 Ophiuchi (data from 
Handler et al, 2005) and right: v Eri (data from Jerzykiewicz et al., 2005) 



again (3q and f3\ are known from the stellar model; 



A> 



K(x) dx ; j3 e 



K(x){x — x c ) dx 



The knowledge of S\ and S2 then yields J?o and fl\ . A rotation gradient in the 
envelope with fl c j Q e = 3.6 is obtained. 

In addition, the seismic modelling of the detected axisymmetric modes favors 
a core overshooting distance of ~ 0.1 pressure scale height (H p ) rather than 
while an overshoot of 0.2 H p is rejected. 



6 Ophiuchi 

is also a main sequence ~ 9 Mq star with an effective temperature T e // ~ 
22900 K. Three multisite campains seismic observations and data analyses reveal 

7 identified frequencies: the radial fundamental t — (pi); one triplet I = l(pl) 
and 3 components (m = —1,1,2) of a quintuplet £ = 2 (gl) (Handler et al, 
2005). A seismic analysis led Briquet et al. (2007) to conclude that the case of 
9 Ophiuchi is similar to HD129929. The detected modes do not provide strong 
constraint about the rotation of the convective core. On the other hand, unlike 
HD129929, the data for 8 Ophiuchi are compatible with a uniform or a quite 
slowly varying rotation of the envelope. The convective core overshoot distance 
is found to be (0.44 ± 0.07) H p This is a much larger amount than found for 
HD129929. Whether this difference must be related to the fact that Ophiuchi 
seems to rotate more than 10 times faster than HD129929 remains an open issue. 
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v Eri 

is a very interesting case as it oscillates with 3 triplets £ = 1 (gi, p\, P2), one 
radial mode p\ and one £ — 2 component. Seismic studies show that the detected 
modes are able to probe the rotation of the core, which is rotating faster than the 
envelope (Pamyatnykh et al., 2004, PHD04; Ausseloos et al., 2004; Dziembowski 
& Pamyatnykh, 2008 (DP08); Suarez ct al. 2009). DP08 further assumed a linear 
gradient as a transition (in the /1 gradient zone) between the uniform fast rotation 
fi = f2 c of the core and the uniform slow rotation of the envelope fl = fi e above 
the /i-gradicnt region. They find a ratio Q c / Q e = 5.3—5.8. Model fitting based on 
the 3 axisymmetric £ — 1 modes yield an extension of the mixed central region of 
0.1-0.28 H p above the convective core radius depending on the adopted chemical 
mixture and metallicity value (DP08; Suarez et al., 2009). 



Table 1. Overshoot versus rotation rate for several stars from seismic analysis. Veq 
is the derived the equatorial velocity, a ov the overshoot parameter, Q inner / Q env the 
ratio of the the rotation rate in the inner layers to that of the surface, Z the metallicity. 
The modellings assume a Grevesse-Noels mixture except for 12Lac. 

/3 Cep Vcq (km/s) a ov ft„„ er /fl e „„ Z ref 

0.1 ± 0.05 n Q .2/Osurf ~ 3.1 0.019 ± 0.003 (I) 
0.44* ±0.07 env. unif. rotation 0.012 ±0.003 (2) 
0.15 ± 0.05 Q c /f2 env ~ 5.5 - 5.8 0.0172 ±0.0013 (3) 

<0.4 n c /n e nv~1.8-5 0.01-0.015 (4), (5) 

*Asplund mixture 

(1) Dupret et al, 2004 (2) Briquet et al, 2007 (3) Pamyathnyck et al., 2004, (4) DP08 

(5) Desmet et al, 2009 



12 Lac 

Several frequencies have been detected for this star (Handler et al., 2006) but 
only 4 of the detected frequencies correspond to identified (£, m) modes (Desmet 
et al, 2009). Only 2 successive components of one £ = 1 triplet are known which 
is not enough to provide information on the inner/surface rotation ratio. One can 
use as an additional information the equatorial surface value, v eq — 49 ± 3 km/s 
as derived by Desmet et al. (2009) . One needs the stellar radius which is derived 
from a seismic modelling of the star. The resulting seismic model and its radius 
depend on the radial orders identified for the modes (DP08 and Desmet et al., 
2009). Second order (centrifugal) effects on the frequencies must also be taken 
into account as the rotation for 12 Lac seems to be fast enough as recognized by 
DP08. Taking then a value for the stellar radius in the broad range R — 7 — 9i?© 
, the equatorial surface value, v eq = 49 ± 3 km/s and the observed splitting of 
1.3032/iHz yields a ratio £2i nner / Q sur f in the range [1.8 — 5] definitely indicating 
a non rigid rotation. There is not yet an agreement concerning the radial order 
of the identified modes but the triplet seems in any case to be of mixed nature 
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and therefore able to probe the core rotation. DP08 did not consider overshoot 
and Desmet et al. (2009) found that core overshoot must be smaller than 0.4 
Hp. 

Summary 

These studies lead to the conclusion that a few rotationally split modes can 
provide important information about internal rotation and core overshoot of 
/3 Cephei stars if the modes are identified, enough precise measurements are 
obtained and the age of the star is such that excited modes have mixed g, 
p nature. Trying to disentangle overshoot and rotation effects on core element 
mixing is only starting with a measure of their relative magnitude as is illustrated 
in Table 1. As emphasized by DP08, in that respect, seismic modelling of fast 
rotators are needed. Once the size of the mixed core and the ratio of core to 
surface rotation are reliably determined, the next issue is to estimate , what 
part in the seismically measured extension of the core, d ov , comes from convective 
eddy overshooting and what part comes from other transport processes such as 
those induced by rotation. 

4.2 Splitting asymmetries : distorsion 

The splitting asymmetry, A m (Eq l35|) . for acoustic modes is mainly due to the 
oblateness of the star caused by the centrifugal force although for low radial 
order modes, the Coriolis contribution is also significant. FigJT] represents the 
normalized splitting asymmetries: 

Tl m = A m /r/ ,i, (48) 

for the I = l,pi and g\ modes in function of the scaled frequency y = vt. n ,o/ ^0,1,0 
where ^0,1,0 is the frequency of the radial fundamental mode. lZ m is plotted for 
9 Ophiuchi, HD129929 and v Eri. The same quantities for 8.2 Mq stellar mod- 
els are also represented. The models have been computed with CESAM2k code 
(Morel, 1997) assuming standard physics (Lebreton et al, 2008; Goupil 2008) 
including a core overshooting distance of 0.1 H p and an initial hydrogen abun- 
dance X = 0.71 and metal abundance Z = 0.014. The evolution of the selected 
models is represented by the central hydrogen content X c from 0.5 to 0.2. The 
frequencies have been computed using a second order perturbation method and 
an adiabatic oscillation code WAR(saw)M(eudon) adapted from the Warsaw's 
oscillation code (Daszyhska-Daszkiewicz et al., 2002). For each model, two sets of 
frequencies are computed assuming a uniform rotation corresponding to v = 30 
km/s and v — 10 km/s respectively. These sequences of models do not represent 
true evolutionary sequences as in realistic conditions, the rotation changes with 
time and can be non uniform. They however illustrate the evolution of the asym- 
metry when a mode changes its nature during evolution, from pure p mode to 
mixed p and g mode for instance. Indeed pure g modes have small asymmetries 
compared with pure p modes because they have much smaller amplitude in the 
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Fig. 7. left: Scaled asymmetries lZ m 10 3 for ! = ln=l (top) and n — — 1 (bottom) 
modes in function of the m = frequency scaled by the radial fundamental mode 
frequency. The open dot (resp. full dot, full square) represents the observed asymmetry 
for 6 Oph, (resp. for V386 Cen, v Eri). The solid (resp. dashed) line corresponds to 
v = 30 km/s (resp. 10 km/s) 8.2Mq models. The central hydrogen content X c is 
decreasing toward the right, top right: Kernels K2(x) for splitting asymmetries of 
£ = 1, n = 1 (pi) mode (solid line) and £ — 1, n — — 1 (gl) mode (dashed line) for 
model with Xc — 0.35. The abscissae is the normalized radius. bottom right Run of 
the normalized Brunt- Vaissala profile N 2 r/g for the corresponding model with r/R. 
from Goupil & Talon, 2008) 



outer envelope where distorsion has its most significant effect. This is illustrated 
in FigJT] In a perturbation description, one finds that lZ m is a second order ef- 
fect proportional to il 2 (DG92; Goupil et al, 2000; Goupil, 2009 and references 
therein) . The variation of TL m with the scaled frequency y (ie with stellar evolu- 
tion) is similar for the v = 30 and v = 10 km/s sequences of models but TZ m is 
roughly 9 times (ie ratio of fi 2 ) larger for v = 30 km/s models than v = 10 km/s 
models. For pure p modes, the asymmetry amounts to TZ m ~ 0.8 10~ 3 whereas 
for pure g modes it almost vanishes. lZ m for the £ = \,n — 1 mode decreases 
for older models (larger y). The reverse happens for the £ = l,n = — 1 mode. 
The reason is that for young models, £ — 1, n — 1 and n = — 1 modes are pure 
p and g modes respectively. When the model is more evolved, these 2 modes 
experience an avoided crossing and exchange their nature. From a perturbative 
approach, one derives: 



(49) 



A m = vi.nfl I ft (x) K 2 (x) dx 



where U 2 = f} 2 /(GM / R?) and x — r/R the radius normalized to the surface 
radius. K2(x) depends on the centrifugal perturbation part of pressure and den- 
sity as well as the differential rotation fi(x) and the mode eigenfunction. Figj7] 
shows K2(x) in function of the normalized radius x = r/R for £ = l,n = 1 (pi) 
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mixing (solid) line and a model including a 0.1 H p overshoot instead (dashed line). 
(from Montalban et al. 2008) 



and I = 1, n = —1 (g\) modes for the v = 30 km/s, 8.2 M© model with X c = 0.5. 
The inner layers contribute to the asymmetry of £ — 1 , g\ multiplet in contrast 
with the £ = l,pi multiplet for which the kernel K 2 is concentrated toward the 
surface layers. The asymmetry of the £ = l,g\ multiplet is sensitive to the inner 
maximum of the Brunt- Vaissala frequency, arising from the /x-gradient, which 
contributes negatively to K^. As the negative contribution is very localized, it 
decreases the asymmetry only slightly compared to a pure p mode for a uniform 
rotation. However, one can expect a larger decrease in case of a rotation faster 
in the inner regions than the surface. 

Theoretical estimates seem to disagree with observed asymmetries deduced 
from £ = 2 modes for Ophiuchi (Briquet et al., 2007) and v Eridani for £ = 
l,p 2 (Dziembowski & Jerzykiewicz, 1999, PHD04, Suarez et al. 2009). Is the 
disagreement real? The question has some relevance as the asymmetry values 
are only marginally above the observation uncertainties. Or can it be that the 
observed frequencies do not belong to the same multiplet as suggested by DP08 
for v Eri? 



4.3 Axisymmetric modes: mixing 

Rotationally induced mixing of chemical elements changes the structure and in 
particular affects the Brunt- Vaissala frequency N at the border of the convective 
core. As a consequence, at a given location in a HR diagram corresponding to an 
observed star, one can find several models with different structures and therefore 
likely different values of the mode frequencies including axisymmetric modes 
which can then be used as diagnostics for mixing. 
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Uniform and constant diffusion coefficient D t : 

Montalban et al. (2008) and Miglio et al. (2008) investigated the effect of tur- 
bulent mixing on a g-mode frequency spectrum and the ability of such modes to 
probe the size of stellar convective cores. They assumed a constant in time and 
uniform in space global diffusion coefficient D t — D e ff + D v in Eq.[T]above. The 
constant value for D t is chosen so as to correspond to the value near the con- 
vective core provided by a Geneva stellar model including rotationally induced 
mixing. This is valid for g-modes which have most of their amplitude there (see 
Miglio et al., 2008) The model is a mid main sequence (X c = 0.3) fO M© with 
D t = 7 10 4 cm 2 /s chosen to correspond to a rotational velocity v= 50 km/s. 

FigH] shows the Brunt- Valssala, frequency (TV) profile for a model with tur- 
bulent chemical element mixing and a model with no turbulent chemical element 
mixing but including instead core overshoot assuming an overshoot distance of 
0.1 H p . Differences can be seen at the edge of the convective core. The Briint- 
Vaissala, frequency of the model with turbulent mixing behaves more smoothly 
in the /i-gradient region above the convective core than for the model computed 
with no turbulent mixing but with an overshoot distance of 0.1 H p . From Geneva 
code calculations, the evolution of the rotation profile leads to a core to envelope 
ratio of 1.6. The differences between the two profiles arising at the edge of the 
convective core cause significant changes on frequencies of g-modes and mixed 
modes. The frequency separations A n j = z^,n,o — ^e,n~i.o differ by a few /iHz 
for radial order n = — 1 and n = — 2 , £ = 2 modes between the model with over- 
shoot 0.1 H p and the model with turbulent mixing (FiglH]). At higher frequencies 
for pure p-modes, no differences in A n j are seen when adding turbulent mixing 
or not. 



Rotationally induced diffusion coefficient 

In this section, we consider stellar models which are computed with the Toulouse- 
Geneva evolutionary code which includes the coupling between rotationally in- 
duced mixing and momentum transport (Eq|T]and Eq.2 above) as described by 
Talon (2008). The rotational evolution of the star begins from solid body when 
the core is still radiative, shortly after the star leaves the Hayashi track. A 8.5 M Q 
mass has been chosen so that the models evolve through the HR diagram to a 
location where the star 9 Ophiuchi is expected (logL/L© = 3.73, T cff = 4.35). 
This corresponds to a mid main sequence model, V15, with a central hydrogen 
content X c — 0.3. The evolution has been initiated with a uniform rotational 
velocity v = 15 km/s on the pms; the rotation profile then evolves to strongly 
differential rotation so that V15 has a surface velocity of v = 48.2 km/s and a ra- 
tio f2 core /f2 sur f — 1.6 when crossing the 9 Ophiuchi location in the HR diagram 
at an age of 19.65 Myr. 

The diffusion coefficient, D t , depends on the meridional circulation velocity 
and the local turbulence strength, ft varies with depth and evolves with time 
as illustrated in Fig|9j The D t profile is represented for 3 models with ages 0.5 
Myr, 1 Myr and 1.5 Myr built assuming an initial 15 km/s velocity on the pms. 
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The rotation evolving from uniform to strongly differential rotation causes a 
relaxation toward a stationary profile which persists with only an ajustement 
due to expansion and contraction with evolution (Goupil & Talon 2008). 
Effect of rotationally induced mixing on the structure is significant at the edge 
of the convective core as emphasized in Fig(S] where we compare the squared 
Brunt- Vaissala profile, N 2 , in the vicinity of the edge of convective core for 
model V15 and a model Vq which includes neither rotationally induced mixing 
nor overshoot. Inclusion of rotationally induced mixing leads to the model V15 
which shows a narrower maximum of Brunt- Vaissala profile at the edge of the 
convective core compared with that of Vo . 

To illustrate the impact of such a difference on the oscillation frequencies, we 
compare low radial order frequencies of the models V15, and Vo- Modes pi, j»2, 
g\ for these models have amplitudes near the edge of the convective core. FigfTTI 
shows that this can result in significant frequency differences for the same mode 
easily detectable with C0R0T observations. The frequencies of these modes are 
quite sensitive to the detail of the Brunt- Vaissala profile in this region. This 
means that some care must be taken when computing these frequencies and 
drawing conclusions. The frequencies of these modes are indeed sensitive not 
only to the physics but unfortunately also to the numerics which can be quite 
inaccurate in this region of the star. 

The sign and magnitude of 8ui = uiy\§ — ^vo ar e dependent on the mode 
when it has amplitude in the regions where the nonrotating model and the 
model with rotationally induced mixing differ. We consider here, as in Sec. 3. 5, 
only the effect of rotationnally induced mixing on the spherically symmetric 
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Fig. 10. Model Vis with a surface rotational velocity v=48.2 km/s. top: Profile of 
of the normalized Brunt- Vaissala, frequency as defined by iV 2 r/g in function of the 
normalized radius r/R. Bottom: rotation profile normalized to its surface value. The 
core to surface ratio for the rotation rate then is 1.6 (from Goupil & Talon 2009 in 
prep.). 



structure. Differences in the structure of the model V15 which includes rotation- 
naly induced mixing and the model Vq which does not result in differences in 
the eigenfrequencies which we note Suj = uiq,o — ^o.i7=o- 

The structure of the models V15 and Vq indicates that pa and its derivative, 
the gravity go, the density pa are not significantly modified compared to the 
derivative of the density. FigfTTIshows that the largest difference S(d log pa/d log r) 
arises near the convective core. Accordingly from Eq l47l one expect larger fre- 
quency differences 5lo for mixed modes compared to p-modes. This is what is 
observed in FigfTTl As explained in Sect. 3. 5, with the help of the integral relation 
for 5uj, the frequency differences for high frequency (i.e. pure) p-mode is small 

and positive. For lower frequency mixed modes, 8 1 " p " ) can be positive and 
the frequency difference can be large and negative as illustrated in FigfTTl 



5 Cubic order versus latitudinal dependence 

It has been known for a long time that latitudinal variations of the rotation 
rate generate departures from linear splitting. On the other hand, a fast uniform 
rotation can generate cubic order corrections to the frequency of non axisym- 
metric modes which also cause departure from linear splitting. The latitudinal 
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Fig. 11. left top: Zoom of Brunt- Vaissala frequency profile in the vicinity of the edge 
of convective core in function of the normalized radius r/R for model V15 (dashed line) 
and model Vb (solid line). The local maximum of N r/g corresponds to a nonzero fi- 
gradient. It decreases more sharply in presence of rotationally induced mixing because 
mixing results in smoothing the ^-gradient. Left bottom: Differences 8v = Suj/(2n) 
between frequencies computed from model Vb (no rotationally induced mixing included) 
and model V15 for £ = (solid line) and £ — 1 (dashed line) m = modes in uilz (from 
Goupil & Talon, 2008). right: dlnp/dlnr in function of the normalized radius r for 
model V15 (dashed line)and model Vb (solid line) right top: from center to surface, 
right bottom: in the central region (from Goupil & Talon 2009, in prep.) 



correction to the linear splitting is proportional to the fl gradient whereas cubic 
order effects, as their name indicate, are proportional to J? 3 . It is expected that 
the dependence of these corrections with the frequency differs when it is due to 
latitudinal differential rotation or to cubic effects. 

Low mass stars are known to be slow rotators. Indeed due to their outer 
convection zone, they undergo magnetic braking. Due again to their outer 
convection zone, obervational evidences exist for surface latitudinal differen- 
tial rotation. Hence for these stars, the averaged rotation rate fl is small and 
Af2 — Q e qua — ^poie, the difference between the rotation rates at the equator 
and the poles, can be large (25% -30% for the Sun, between 1% and 45% for a 
star like Procyon, Bonanno et al. 2007). One therefore expects that latitudinal 
corrections to the splittings dominate over cubic order ones which are negligi- 
ble. On the other hand, more massive stars on the main sequence have shallower 
convection zones which even disappear above ^3-5 M®. These stars usually are 
fast rotators with a radiative envelope which may or may not be in latitudinal 
differential rotation. For these fast rotators, one can wonder what is the minimal 
latitudinal shear which dominates over cubic order effects and can therefore be 
detectable. Here we quantify this issue with the help of a polytropic model with 
index 3. The constants characterizing the polytrope are taken to correspond to 
model A considered in Sect l4.3l We establish first the splitting correction due to 
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latitudinal differential rotation. This is then compared with the splitting correc- 
tion arising from cubic order effects as derived by previous works. We assume a 
rotation velocity of 100 km/s. 



5.1 Latitudinal dependence 

Hansen et al. (1977) derived the expression for the rotational splitting of adia- 
batic nonradial oscillations for slow differential (steady, axially symmetric) rota- 
tion J?(r, 9) and applied it to numerical models of white dwarfs and of massive 
main sequence stars assuming a cylindrically symmetric rotation law. In the so- 
lar case, the effects of latitudinal differential rotation on theoretical frequencies 
were investigated by Gough & Thompson (1990), Dziembowski & Goode (1991, 
DG91) and DG92 who also considered the case of S Scuti stars. 

In order to be able to compute the splittings from EqfT7land Eq|22l one must 
specify a rotation law. It is convenient to assume a rotation of the type: 



Q{r,6)= ]T n 2s {r) {cost 



\-s 



(50) 



s=0 



where 9 is the colatitude and we take s ma x = 2. The surface rotation at the 
equator is D(r = R,9 = n/2) = fio(r = R). Note that in the solar case, f2 2 , fl^ 
are negative and the equator rotates faster than the poles (DG91, Schou et al, 
1998). As shown in Appendix, inserting Eq[50] into Eq|22l yields the following 
expression for the generalized splitting (Eg 1861 in Appendix): 



s=1 



S m = f (2 (r) K{r) dr + V m 2s H S {Q) 
with K(r) defined in Eql25l and 



(51) 



H a {Q) 



1 /■■" 
- j J fl (r) [Rs (C ~ 2fr& + H{A - 1)) + Q s e h )]por 2 dr (52) 

where R s and Q s depend on Q 2 , ^4 and A — £(£ + 1) and are given by Eg J5H 
and Eg 1861 (Appendix) respectively. 

a) Uniform rotation In that case, i?(r, 9) = Qq, Vr, 9; i?2, ^4 = i.e. Rj, Qj = 
for j = 0, 2 hence H m .j — 0. One recovers the well known expression: 



S m — ^0 ft 

where, for later purpose, we have defined 

/3= I K{r) dr 
Jo 

r-fl 



(53) 



£ - 2£ r & + (A - 1) £ Por 2 dr 



(54) 
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This is usually rewritten as : 

S m = f2o (Cl — 1) 
where Cl is the Ledoux (1951) constant 



C L 



1 



2£r& + ti 



por dr 



b) Shellular rotation then f2(r 7 9) — £2o(r) and s r 
f2 2 , Mi = ie i?j = Qj : = for j = 0, 2 and 



s -- 1 - 



flo(r) 



£-2&.&+j1$ Por 2 dr 



0; again here: 



(55) 



c) Latitudinally differential rotation only In that case, Qij-,3 = 0,2 are 

depth independent and R s and Q s are constant and 



s=2 



S m = Q P + n Y, ™ 2s (R S {Q) + Q S {Q) 7) 

s=0 

with j3 defined in EqJSH and 

1 f R 

J Zh Por 2 dr 

1 Jo 



(56) 



7 = 



For a triplet £ = 1, m = 1 (-4 = 2) then 

Si = Q (3 + n {R{Q) p + Q[Q) 7 ) 
with (using Eq|5H and Eq[ 



R{ f2)=J2R s (n)= 1 ~%L + |§ 

^ 5 J 2 7 J / 

Q(r2 ) = gQ s( r2) = -|^ 



(57) 

(58) 
(59) 



In the solar case, j3 ~ — 1 and |/3| >> |7| for the excited high frequency 
p-modes. 



Si « -/?,, (1 



1 r? 2 3 /? 4 

5 I 2q 7 J/q 



(60) 



With f2 2 /f2a = -0.127, r2 4 /^o = -0.159 (from DG89), one obtains a departure 
from linear splitting of \Si/Oq + 1\ = 0.093 i.e. a w 10% change in the solar case. 
For upper main sequence stars, excited modes are around the fundamental radial 
mode and may be mixed modes with |/3| ~ I7I ~ 1/2. This leads for instance to 
\Si/{2q + 1/2 1 « 5% for i? 2 /^o an d fi^/Qo equal to 1/5 of the solar values. 



/ 1 AH\ 
Si(lat) = f2 l3(l + ~ — J (61) 
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5.2 Latitudinal dependence versus cubic order effects 

Let assume on one side a pulsating star uniformaly rotating with a rate i?o 
high enough that cubic order (0(J?q)) contributions are significant. On the other 
side, one also considers a model rotating with a latitudinally differential rotation 
(uniform in radius). One issue then is which one of these two effects dominate 
over the other one since the cubic one is 0(f2 3 ) whereas the other one is 0{AQ) 
? For stars other than the Sun, one can simply assume the rotational latitudinal 
shear Ail — il 2 with J? 4 = and fi[6) = flo + Afl cos 2 6. For £ = I modes, 
Eq|5"Tl becomes 

1 AG> 

Expressions for the frequency correction (in rad/s) for cubic order effects 
assuming a uniform rotation has been derived by Soufi et al. (1998). Part of 
the cubic order effet is included in the eigenfrequency ujq^o and therefore is also 
included in second order coefficients which indeed involve u>o,Q. Another part of 
the cubic order effects is included as an additive correction to the frequency. 

Frequency up to 3rd order were computed for models of 5 Scuti stars by 
Goupil et al (2001), Goupil & Talon (2002), Pamyatnykh (2003), Goupil et al. 
(2004). Karami (2008) rederived the cubic order effects following Soufi et al.'s 
approach and Karami (2008, 2009) applied it to a ZAMS model of a 12 M Q 
/? Cephei star. He found that cubic order effects are of the order of 0.01% for 
a I — 2, n — 2 and 0.5% for a n=14 mode for a 100 km/s rotational velocity. 
Values of the third order additive correction to the frequency were listed for 
£ = 1 p- modes of a polytrope of index 3 by Goupil (2009). 

Here we write the splitting under the form: 

S m (cubic) = On (3 + f2 ( — J T\ m \ (62) 

where the last term represents the full cubic order contribution with Co is the 
normalized frequency of the nonrotating polytrope and SIq = Qq/Qk- 

Tabf2] lists the value of the dimensionless coefficients T\/a\] and — /3, —7 for 
£ — 1 modes for a polytrope with a polytropic index 3. The coefficient Ti/<7 2 
remains nearly constant with increasing frequency for frequencies above co > 10 
i.e. for p modes For 00 > 10 (p-modes), — (3 ~ 1 and Ti/cr 2 « —0.09. The splitting 
is decreased by a latitudinal dependence with AQ < whereas it is increased 
by cubic order effects Ti/(3 > . In absolute values, the effect of latitudinal 
differential rotation on the splittings then dominates over cubic order effects 
whenever: 

1^1 > 0.45 nl 

J/o 

For model A and a rotational velocity 100 km/s, fi = 0.174 then \^f-\ > 1.36% 
For a faster rotator with for instance 200 km/s, the latitutinal shear must be 
larger i.e. \jf-\ > 5.45%. 
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For the slowly rotating f3 Cep stars considered in Sect. 4 above (v < 50 km/s), 
cubic order effects in the splittings can be neglected in front of latitudinal effects 
equal or larger than 0.34%. At this low level, both effects are comparable to the 
observational uncertainties (0.1%). 



Table 2. Coefficients assuming a uniform rotation for a polytrope with polytropic 
index 3 and adiabatic index 7 = 5/3. The squared frequency <tq is the dimensionless 
squared frequency u) /(GM/R ). Spherical centrifugal distorsion of the polytrope has 
not been included. 

1=1 
n ^ UTo Ti T 2 Y 1 Y 2 Ti/<7g -P ~ 



-7 


0.22 


0.479 


0.417 


0.008 


0.012 


-0.018 


0.592 


0.521 


0.456 


-6 


0.28 


0.476 


0.419 


0.005 


0.015 


-0.023 


0.462 


0.524 


0.450 


-5 


0.37 


0.473 


0.422 


0.001 


0.020 


-0.029 


0.351 


0.527 


0.441 


-4 


0.52 


0.469 


0.425 


-0.004 


0.026 


-0.039 


0.254 


0.531 


0.431 


-3 


0.78 


0.466 


0.427 


-0.013 


0.038 


-0.056 


0.164 


0.534 


0.410 


-2 


1.28 


0.466 


0.428 


-0.024 


0.059 


-0.089 


0.073 


0.535 


0.386 


-1 


2.51 


0.473 


0.422 


-0.035 


0.106 


-0.159 


-0.025 


0.528 


0.269 


1 


11.37 


0.029 


0.777 


0.877 


2.890 


-4.335 


0.024 


0.970 


0.025 


2 


21.49 


0.034 


0.773 


0.864 


5.802 


-8.703 


-0.034 


0.966 


0.028 


3 


34.83 


0.033 


0.773 


0.851 


9.624 


-14.436 


-0.063 


0.966 


0.027 


4 


51.39 


0.031 


0.776 


0.840 


14.340 


-21.511 


-0.077 


0.969 


0.026 


5 


71.15 


0.027 


0.778 


0.832 


19.940 


-29.909 


-0.084 


0.973 


0.025 


6 


94.09 


0.024 


0.781 


0.826 


26.414 


-39.621 


-0.088 


0.976 


0.023 


7 


120.19 


0.021 


0.783 


0.821 


33.757 


-50.635 


-0.089 


0.979 


0.022 


8 


149.43 


0.019 


0.785 


0.817 


41.964 


-62.946 


-0.089 


0.981 


0.020 


9 


181.81 


0.017 


0.787 


0.814 


51.032 


-76.548 


-0.089 


0.984 


0.019 


10 


217.32 


0.015 


0.788 


0.811 


60.958 


-91.437 


-0.089 


0.985 


0.018 


11 


255.94 


0.013 


0.789 


0.809 


71.739 


-107.609 


-0.088 


0.987 


0.017 


12 


297.67 


0.012 


0.790 


0.807 


83.375 


-125.062 


-0.087 


0.988 


0.017 


13 


342.51 


0.011 


0.791 


0.805 


95.862 


-143.793 


-0.087 


0.989 


0.016 


14 


390.44 


0.010 


0.792 


0.804 


109.201 


-163.802 


-0.086 


0.990 


0.015 


15 


441.47 


0.009 


0.793 


0.803 


123.392 


-185.087 


-0.085 


0.991 


0.014 


16 


495.59 


0.008 


0.793 


0.802 


138.432 


-207.648 


-0.085 


0.992 


0.014 


17 


552.80 


0.008 


0.794 


0.801 


154.323 


-231.484 


-0.084 


0.993 


0.013 


18 


613.09 


0.007 


0.794 


0.800 


171.064 


-256.595 


-0.084 


0.993 


0.013 


19 


676.47 


0.006 


0.795 


0.799 


188.655 


-282.982 


-0.083 


0.994 


0.012 


20 


742.93 


0.006 


0.795 


0.798 


207.097 


-310.645 


-0.083 


0.994 


0.012 


21 


812.46 


0.006 


0.796 


0.798 


226.389 


-339.584 


-0.082 


0.995 


0.012 


22 


885.08 


0.005 


0.796 


0.797 


246.534 


-369.800 


-0.082 


0.995 


0.011 


23 


960.78 


0.005 


0.796 


0.797 


267.530 


-401.295 


-0.082 


0.995 


0.011 
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6 Conclusions 

We have seen along this review that several efficient seismic tools can be de- 
signed to obtain valuable information on the internal structure and dynamics of 
main sequence massive stars which oscillate with a few identified modes. Iden- 
tification of the detected modes requires a high signal to noise which is made 
available due to the large amplitudes of these opacity-driven modes. On the other 
hand, these stars oscillate with low frequencies lying near/in the dense part of 
the spectrum where p modes, mixed modes and g modes can be encountered. 
While this is a great advantage in order to probe the inner layers of the star, 
resolution and precise measurement of quite close frequencies in a Fourier spec- 
trum requires very long time series. This explains the yet still small number of 
f3 Cephei stars for which a successful seismic analysis has been obtained, de- 
spite the appealing prospects that a better knowledge of their structure bring 
up valuable constrains on their still poorly understood life end as supernovae. It 
is expected that the space experiments CoRoT (Michel et al. 2008) and Kepler 
(Christensen-Dalsgaard et al, 2008) will increase the number of O-B stars for 
which fruitful seismic analyses can be carried out as well as possibly enlarge the 
sample to fast rotators. Mode identification can be at first difficult to perform 
for fast rotators but some of these fast rotating stars might also show oscilla- 
tions of solar-like type which characteristics could help the mode identification. 
This interesting perspective has recently emerged with the discovery of the first 
chimera star with the CoRoT mission (Bclkacem et al, 2009). 
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Appendix : Differential rotation 

The expression for the mode splitting of adiabatic nonradial oscillations due to 
a differential rotation fi{r, 9) can be put into the compact form (Hansen et al. 
(1977), DG91, DG92, Schou ct al. (1994a,b), Pijpers (1997), CD03): 

6oj m =m / K m {r,6) f2(r,6) dQdr (63) 



Jo Jo 
where /C™ is called rotational kernel: 



k^) = -^^ M )x 



J^[(\&\ 2 -(CM + cc))\Yr\ 2 (64) 
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where the spherical harmonics Y™ are normalized such that 

(Fr')*(M) Y e m (e,4>) ^ - 6 t ,t>6 m , m . 

where dfl = sin 9d9d<p is the solid angle elemental variation and Se^' is the 
Kroenecker symbol. Mode inertia I is given by 



(\&\ 2 + A\£ h \ 2 ) P0 r 2 dr (65) 



with A = £(£+l). 

It is convenient to assume a rotation of the type: 



:?.", 



fl(r,0)= £ r? 2s (r) (cos( 



s2s 



(66) 



where is the colatitude. Eg |6"51 becomes: 



SiOr, 



(l^| 2 ~ (£& + cc)) 5 S + |a| 2 (Bi + S 2 ) 



(67) 



por 2 dr 



where we have defined 



S s 



with /^ = cos 9 and 



y/T (cos i 



\^>' 



dQ 

47T 



^ 2s |r/"(^,0)r^ 



(68) 

(69) 
(70) 



The term in |£/j| 2 requires a little care. Consider first Bl. Integration by part 
leads to 



Bi = J (VbYT ■ V H Y e m )(co S 9) 2s ^ 



B\ 



dQ 

An 



17 



(71) 



(V?^" 1 ) (co S e) 2s + (V H Y e m ) • V ff ((cos0) 2s ) 



Recalling that V^F/™ = -A Y e m , one gets 



Bl = A S„ 



1 



9F™ d(cos0) 2s ) 



47T 



96* d6» 

where cc means complexe conjugate. Again an integration by part yields 

d(cos(9) 2s )i dO 



bi - a s. + UprrZ 



sine 



dd 



One finally obtains 



Bl = A S„ 



(72) 



(73) 



(74) 



(2s - l)S.-i - (2s + 1)S S 

Turning to the second term B 2 in Eq |68l an integration by part yields 

B 2 = -(2s+l) S s (75) 
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Inserting expressions EqJTH and [73] into Eq. [551 one obtains 



SUm = 



E 

s=0 



^2s{r) K„ hS {r)dr 



with 



and 



K m ,.(r) = K (r) S s -j Po r 2 \^ h \ 2 s [{2s - l)5._i - (2s + 3)5 S 



K{r) = -j [\ZA 2 {CAh + ce) + \S h \ 2 (A 1) 



PO T 



(76) 



(77) 



(78) 



Expression EqJTT] is equivalent to Eq.25 in DG92. For any s, S s is given by 
a recurrent relation (Eq.31 in DG92). Note that 6u m = 5ui- m . Let define the 
generalized splitting 



*Jm. — 



8u m — 8uj- 



5uj„ 



2m 2m m 

We limit the expression for the rotation to s max — 2 i.e.: 

Q(r, 6) = n (r) + Q 2 {r) cos 2 6 + Q A (r) cos 4 6 

then for adiabatic oscillations {£r{r) and £fc(r) are real): 

K mfi (r) = K(r) 



(79) 



1 



K mA (r) = K{r) S x - j & 
K m ,2(r) = K(r) S 2 - \ £ 2 2 



l-55i 



Por 



(80) 



35i - 75 2 



Por 



where we have used <S_i = 0; <So = 1- 

We obtain a formulation for the generalized splittings with a m dependence 
of the form: 

S m = I (fl + Q 2 S l + Q 4 S 2 ^j K{r) dr 
-j j (n 2 {l - 5Si) + Qa 2(35! - 7S 2 )) & p r 2 dr (81) 



One needs Si and S 2 (computed from Eq.31 in DG92): 



Si 

s 2 



1 



4/1-3 
1 3 



-(-2m i + 2A-l) 



2A-1 
4/1-3 



4/1-3 



[ l Si(-2m 2 + 2/l-9) + l] 



4/1-15 2 
The first term in brackets in Eq|HT] becomes 
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For the second term in Eal81| one has: 

f? 2 (l - 55i) + Q 4 2(35i - 7S 2 = f2 (<2o + m 2 Q 2 + m 4 Q 2 ) 
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Collecting terms from Eq|82land Eq|84l the generalized splitting Eq|8T1 takes 
the expression: 



S m = / flb(r) K(r) dr + Y] m 2s H S (Q) 
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For a depth independent rotation law, f2{9), f^2j,j = 0, 2 are depth indepen- 
dent and R s and Q s are constant, then for a triplet I = 1 (/l = 2): 



Si = /?„ /3 + fl (I] fl.(fl)) £ + £ Qs{0)) 
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